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Abstract 

A relaxation process, with the associated phenomenology of sound attenua- 
tion and sound velocity dispersion, is found in a simulated harmonic Lennard- 
Jones glass. We propose to identify this process with the so called microscopic 
(or instantaneous) relaxation process observed in real glasses and supercooled 
liquids. A model based on the memory function approach accounts for the 
observation, and allows to relate to each others: 1) the characteristic time and 
strength of this process, 2) the low frequency limit of the dynamic structure 
factor of the glass, and 3) the high frequency sound attenuation coefficient, 
with its observed quadratic dependence on the momentum transfer. 

PACS numbers:63.50.+x, 61.43.-j, 61.43.Fs 
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In a recent paper W. Gotze and M.R. Mayr |l| adapted the Mode Coupling Theory 
to glassy phase. This theory, as shown numerically for a hard-spheres model, accounts for 
many of the features found in the dynamic structure factor, S{Q,uj), of glasses. Among 
these: i) The existence of propagating excitations (sound waves), with an almost linear 
momentum transfer (Q) dependence of their excitation energy fl{Q), up to Q values that 
are a significant fraction of Qo, the maximum of the static structure factor S{Q). ii) The 
quadratic dependence on Q of the excitations broadening r(Q), and therefore of sound 
attenuation, in) The temperature insensitivity of T{Q). iv) The development in the S{Q, u) 
at large Q ((5/<5o~0.3) of a secondary excitation band, at frequencies below the Brillouin 
peak. The theory also predicts two features not yet experimentally detected: a) A positive 
dispersion of the sound velocity and 6j an intensity "gap" in the low frequency region of the 
S{Q,uj). In spite of the success of this approach, it is still of great interest to investigate 
the physical origin of these phenomena, and, in particular, whether they are related to the 
topological disorder and/or to the anharmonicity of the interatomic potential. 

In this Letter we report a Molecular Dynamics (MD) simulation study of the Q depen- 
dence of the sound velocity in a model monatomic Lennard- Jones glass in the harmonic 
approximation. We show that even in a harmonic glass, by increasing Q there is a positive 
dispersion of the sound velocity, thus proving one of the prediction of the MCT theory [|l|, 
and relating this phenomenon to the topological disorder. Since this dispersion is similar to 
that found in presence of a relaxation process, we attempt to apply a generalized Langevin 
equation with an effective memory function approach to describe the density fluctuations 
dynamics. This formalism allows to account for the ubiquitous Q^-dependence of the high 
frequency sound absorption observed in many glasses by experimenys and by MD sim- 
ulations 0,3. These results suggest to identify this process with the one referred to as 
microscopic or instantaneous relaxation in real systems. 

The investigated systems consist of A^=2048, 10976 and 32000 argon atoms interacting 
via a (6-12) Lennard- Jones potential (e=125.2 K, (T=3.405 A). A standard microcanonical 
MD simulation, performed at decreasing temperatures in the normal liquid phase, is fol- 
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lowed by a fast quench (T^ 10^^ K/s) of the shghtly supercooled liquid down to ~5 K 0. 
Starting from the glass configuration at T=5 K, the atomic trajectories ri{t) (here i=l...N 
is the particle label) are followed and stored for subsequent analysis, and the "inherent" 
configuration {xi}j=i..Ar ||^ at T=0 K is calculated by the steepest descent method. Two 
different procedures have been used to derive the S{Q,uj): (i) The trajectories calculated 
from MD in the glass are used to compute the corresponding time correlation function and 
(ii) The Normal Mode Analysis (NMA) is applied to the inherent configuration. This last 
procedure works in the harmonic approximation, which is obtained by retaining only the 
quadratic term of the interaction potential. From the calculation of the dynamical matrix D, 
one computes the eigenvalues, ujp {p=1..3N is the mode index) and the eigenvectors {ep{i)). 
For large size samples (where a direct diagonalization of D is not feasible) the S{Q,uj) is 
obtained by the method of moments 

In the MD runs, the dynamics structure factor is calculated from the numerical time- 
Fourier transform of the intermediate scattering function F{Q,t), defined as: 

F{Q,t) = l/Ar(Si,exp(zQ-f,(t))exp(-^Q-f,(0))). (1) 

In the harmonic framework and in the classical limit, the one excitation approximation of 
the dynamic structure factor S^'^\Q,uj) is obtained by the normal mode expansion of the 
atomic displacements: 

fiit) =Xi + ^KBT/MT.pep{T)Ap{t)/^p. (2) 

Ap{t) is the amplitude of the p-th normal mode and is characterized by {\Ap{t)\'^) = 1. This 
gives: 

S^'Xq, u) = {KbTQ^/Muj^) T.pEp{Q)5{uj - Up), (3) 

where we have introduced the spatial power spectrum of the longitudinal component of the 
eigenvectors, Ep{Q): 

Ep{Q) = |Ei(Q-ep(0)exp(zg-x,)|l (4) 
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Here, Q=Q/\Q\ and the Debye- Waller factor has been neglected. Differently from the sec- 
ond uj'^S{Q,uj)duj= J iu'^S^^\Q,uj)dij=KBTQ'^ /M) and higher moments sum rules, the 
zeroth moment sum rule for S{Q,ij), J S (Q , ij)duj=S (Q) , does not hold for S^'^\Q,iu), as in 
this function the elastic intensity is missing. Rather, it is useful to define the "inelastic" 
contribution to S{Q): S-{Q)=J S^^\Q,uj)duo. 

Selected examples of the S{Q^uj) calculated for different size systems, and with different 
methods are reported in Fig. 1. As there is a trivial dependence of the S'^^XQ^uj) on T, here 
we report this quantity multiplied by the factor (M/KbT). Figure la shows the S{Q,uj) in 
the low Q range, as calculated from the MD runs for the A^=32000 particles system, while 
Fig. lb shows the intermediate Q range, in the case of the harmonic approximation for the 
A^=2048 system. As a check of consistency, the inset of Fig. lb shows the Q dependence of the 
second moment of {M/ KbT)S{Q, uj) (+) which, according to the sum rule, should be (full 
line). The inset of Fig. la shows, at two Q values, the comparison of the S{Q, uj) (A^=2048) 
calculated either from the MD runs (full line) or in the harmonic approximation (o). The 
S{Q,uj), as derived with the two different methods, are equivalent in the whole considered 
Q range. This indicates that the Newtonian dynamics at T=5 K is truly harmonic, and that 
the results obtained by MD at this T and NMA can be interchanged among each other. 

The general features of the S{Q,uj) reported in Fig. 1 are: i) A Brillouin peak, dispers- 
ing and becoming broader with increasing Q, dominates the spectrum up to Q~10 nm~^ 
(Q/Qo~0.45). ii) At larger Q values a second peak, observed at frequencies below the Bril- 
louin peak, starts to dominate the data. This secondary peak has been already detected in 
many systems: 1) in liquid water, both experimentally and by MD where it has been 
interpreted as a signature of the transverse dynamics; 2) in vitreous silica, by MD, where it 
has been interpreted either in terms of a transverse dynamics or as an evidence of the 
Boson peak |]ll|; and 3) in a hard sphere glass, where it has been theoretically predicted 
|l[]. In this respect, it is worth to underline the striking similarity between the S{Q,uj) 
reported here and those of Fig. 5 of Ref. [jl|. Finally, in), the S{Q,uj) up to Q~5 nm"^ 
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shows a nearly constant and Q-independent intensity below the Brillouin peak frequency. 
This plateau, whose value is \im^_^Q{M/ KbT)S{Q,u)= 74o~0.03 10"^*^ s^/m^, can be iden- 
tified with the plateau observed in the constant uj cuts of the S{Q,uj) at Q larger than the 
Brillouin peak in simulated [0, calculated (Fig 18 in Ref. and measured [|13| systems. 

Figure 2 shows the position of the maxima {Q^{Q)) of the current spectra Cl{Q,u!) 
{=uj'^ /Q'^S{Q,uj)) as a function of Q in the investigated Q range. The inset of Fig. 2 shows 
the peak position VL{Q) and the broadening T{Q) of the S{Q,uj) in the low Q region, as 
measured directly from the spectra of Fig. la. The Q dependence of these parameters 
agrees with the behavior observed in all the glasses investigated so far [||: a linear (square) 
dependece of VL{Q) {V{Q)). The behavior of ^^^Q) closely resembles that of the acoustic 
phonon branches in crystals, i. e. the almost linear behavior in the small Q region, a 
maximum around Qol'^-, and a minimum around Qo- Most importantly one can clearly 
detect a positive dispersion of the sound velocity, as highlighted by the dashed line in Fig. 2. 
This dispersion is better seen in Fig. 3, where the apparent sound velocity, v{ff)=^^{ff)IQ^ 
is reported. The quantity v{ff) (full dots) undergoes a transition from the "low frequency" 
sound velocity Voiff) towards the infinite frequency sound velocity t'oo(<5)- Here Voiff) 



is calculated as Vo{Q) = \J KbT / MSi{Q) (dashed line), and foo(Q) is calculated either as the 
fourth moment of the calculated S{Q,uj) (open points) or from the expression of the fourth 
moment in terms of both the pair correlation function and the interaction potential (full 
line) |16[. The velocity dispersion is observed up to Q~5 nm~^, i. e. i7p~35 cm~^, and an 
opposite dispersion is observed in the region approaching Qo. These observations recall of a 
typicall relaxation scenario. When a relaxation process with characteristic time r is active 
in the system, the transition from Vo to takes place when the condition ujt=1 is fulfilled. 
In the present case, considering that the first transition is at at Qk,2 nmT^ where fi^^lS 
cm~^, the value of r results to be around 0.3 ps. The second and third transitions between 
Vo and foo are observed just below and above Qo as a consequence of the slowing down of 
the dynamics (deGennes narrowing) around Qo- Consequently the whole behavior of v{Q), 
as reported in Fig. 3, can be qualitatively understood in terms of a relaxation process with 



a characteristic time of ps. 

The identification of a such a relaxation process in an harmonic system, suggests to 
use the formahsm that describes the density correlators (j){Q,t)=F{Q,t)/ S{Q) through its 



generahzed Langevin equation |16 



0(Q, t) + ul<P{Q, t) + /*m(Q, t - t')^iQ, t')dt = 



(5) 



where ujI=KbTQ'^ /MS{Q) and m(Q,t) is the " memory function" . This equation has been 
rigorously derived for ergodic systems, but, as recently shown in the framework of the MCT 
lU, it can be still applied in the non-ergodic glassy phase making the the following substi- 
tutions: 0(g,t) ^ <P\Q,t) = {<P{Q,t) - - /J, ul{Q) ^ ^ujI^^=KbTQ^ / MS^iQ) 

and m{Q,t) mfj,{Q,t)=m{Q,t) —u^f^/ {1 — f^). Therefore, the "vibrational" dynamics 
of interest is now described by the correlators (f)'{Q,t), which is obtained subtracting from 
(j){Q,t) the long time plateau level, whose value is the non-ergodicity parameter f^. It is 
worth to note the subtraction of the constant term S{Q)fg from F(Q,t) is equivalent to ne- 
glect the elastic contribution in the S{Q,u!), and that F'{Q,t)=Si{Q)(f>'{Q,t) is the Fourier 
transform of the S^^\Q,uj) as defined in Eq. (^. The whole dynamic behavior of (j)'{Q,t) is 
now contained in the microscopic contribution to the memory function mn{Q, t). In the case 
of a harmonic system, the function m^{Q, t) can be explicitly calculated from the eigenvalues 
and eigenvectors of the system. Indeed, the Laplace transform (indicated by hats) of Eq. @ 
(after the previously indicated substitutions) and a straightforward algebra gives: 



l-s0'(Q,s) 



(6) 



Then, from an inverse Fourier and a subsequent Laplace transform of Eq. (^, it is easy to 
get an explicit expression for 0'(Q, s) to be inserted in Eq. (|^). This gives: 



p J 



(7) 



This equation provides an explicit expression of the memory function in terms of the system 
eigenstates. Considering that m^(Q,t) is mainly characterized by parameters as its initial 
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value and total area and therefore by a decaying time-scale Tq^Fq/A^ [0, one can 
show that m{Q, s^oo)=A^/s and m{Q, s— >0)=Fq |T6[. Inserting these limiting values into 
Eq. (|^ one obtains: 

-1 



A^ 

Q 



-2 



P 9 



lim 



(8) 
(9) 



It is now easy to identify, through the explicit expression of the zeroth and fourth moments 
of Eq. (I) that EpEp{Q)ujl=vlQ^ and {EpEp{Q)uj;^)~^=vlQ^ confirming that ^l = {vl- 
vDQ'^. The determination of F^ which, being the area of the memory function, coincides 
with the Brillouin broadening in the Ci;r<<l limit is slightly more involved. Using the 
representation 5{x) = l/vr lims_»o + x^), in Eq. (j^) one sees that: 



and comparing of Eqs. (|10]) and (^) one gets : 



(10) 



^-^|^lim5(Q,.;). 



(11) 



Similarly: 



ttM 



KnT 



lim S{Q, uj) 



(12) 



Considering that lim^^o '5'(Q, c^;) is Q-independent (see Fig. 1), these expressions give ac- 
count of the observation in the low Q region -where the Q dependence of Vo and v^o 
is negligible- that the Brillouin peak broadening is proportional to and the r is Q- 
independent. As a check of internal consistency, using the value of Ao previously reported, 
one finds r?a0.6 ps and Fp[cm~^]= 1.35(5^[nm^^]. The value of overestimates the one 
deduced from Fig. 3, and this probably due to the rough expression of r as F^/A^^. On the 
contrary, as shown in the inset of Fig. 2, F^^ from Eq. ([Tl|) is in excellent agreement with 
the one directly derived from the width of the S{Q,uj). 

It is worth to discuss the microscopic origin of the observed relaxation process. A relax- 
ation process can be pictured as the macroscopic manifestation of microscopic phenomena 



associated with the existence of channels by which the energy stored in a given "mode" 
relaxes towards other degrees of freedom. The S{Q*,uj), through the fluctuation-dissipation 
relation, reflects the time evolution of the energy initially stored {t=to) in a Plane Wave 
(PW) of wavelength iTxjQ* . As the PW is not an eigenstate of the disordered system, at 
t > to there will be a transfer of amplitude from this PW towards other PWs of different 
Q values. This process is controlled by the difference bewteen the considered PW and the 
normal modes of the topologically disordered glassy structure. This energy flow takes place 
on the time scale r as derived from Eq. ([12|) , and gives rise to the observed relaxation process 
phenomenology. Consequently, one can speculate that this process is the instantaneous or 
microscopic process empirically introduced to explain the S{Q, u) measured in real glasses 
and liquids by Brillouin light and x-ray Scattering |[18|JT9| . 



We thanks W. Gotze for helpful discussions and critical readings of the manuscript. 
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FIGURES 

FIG. 1. The quantity (M/KbT)S{Q,uj) is reported in absolute value at the indicated Q values for 
a) the very low Q region (iV=32000, MD calculation) and b) the intermediate Q region (A^=2048, NMA 
calculation). The inset of a) shows the comparison (iV=2048) of the MD S{Q,u;) (fuU line) and the NMA 
s'' \q,lu) (o). The inset of b) shows the Q dependence of Af^^', the second moment of {M/KbT)S{Q,lu) 
(+) and its theoretical prediction m[^—Q'^ (full line). 

FIG. 2. Frequency position of the maxima of the longitudinal current spectra {^m{Q)) for different 
indicated size samples. The inset shows in log-log scale the maxima {Vl{Q)) and the broadening {V{Q)) of 
the Brillouin peaks in the low Q region. The full line is the prediction according to Eq. (|l^) and the dashed 
line is the best fit to the data. 

FIG. 3. Sound velocity in glassy LJ Argon: apparent sound velocity (w, full symbols), zero frequency 
sound velocity (uq, dashed line) and ininite frequency sound velocity (wqo, full line and open symbols). 
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